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Abstract 



We consider for the Helmholtz equation the inverse problem of identifying locations and 

certain properties of the shapes of small dielectric inhomogeneities in a homogeneous background 

QQ ' medium from boundary measurements on part of the boundary. Using as weights particular 

background solutions constructed by solving a minimization problem we develop an asymptotic 

(variational) method based on appropriate averaging of the partial boundary measurements. 
^ ' 

. Identification de petites inhomogeneites dielectriques a partir de mesures sur une 

partie du bord 

Resume 

Nous considerons le probleme d'identification de petites inhomogeneites dielectriques a partir 
de mesures (incompletes) sur uniquement une partie du bord. Grace a un resultat de densite, 
nous construisons une fonction dont la transformee de Fourier inverse permet de localiser les 
petites inhomogeneites. 



Version frangaise abregee 



■ Soit fl un ouvert borne de R'', d > 2, de classe C^. Supposons que Q contient m inho- 

k>( I mogeneites,{zj + aBj}JLi, on a est un petit parametre, Bj C R'^ est un ouvert borne et les 

j_j ■ points {zjYj^i verifient les hypotheses (|I|). Solent la permeabilite magnetique fXa et la permit- 

I tivite electrique de forme (||). 

Cette Note concerne le probleme de reconstruction des points {zj}™^^ et des tenseurs de 
polarization Mj des domaines Bj, definis par (^, a partir de la mesure de la derivee normale du 
champ electrique Ea, solution de I'equation de Helmholtz (||), sur une partie Ti CC dfl. Grace 



a un resultat de densite, etabli dans la Proposition 2.1, et a la formule asymptotique (^, nous 
reduisons ce probleme inverse au calcul de la transformee de Fourier inverse de la fonction A^, 
definie dans (||). La fonction test Wa, utilisee dans (^), pent etre construite numeriquement en 
resolvant le probleme de minimisation (^. Cette fonction verifie les proprietes d'approximation 
ennoncees dans le Lemme 2.1. Cette Note generalise la methode introduite dans ||l| aux situa- 



tions oil on ne disposerait pas de la mesure de la derivee normale du champ electrique Ea sur 
tout le bord dil. 
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1 Introduction 

Let n be a bounded C^-domain in R'', d> 2 and v be the outward unit normal to 90. Assume 
that contains a finite number of inhomogeneities, each of the form Zj + aBj, where Bj C R'^ 
is a bounded, smooth domain containing the origin. The total collection of inhomogeneities is 
= ^JLi{zj + ceBj). The points Zj € — 1,. . . ,m, which determine the location of the 
inhomogeneities, are assumed to satisfy the following inequalities: 

\zj - zi\> Co > 0,y j I and dist(zj,(917) > Co > 0,V j. (1) 

2d\n\ 

This assumption implies m < Assume that a > 0, the common order of magnitude of the 

diameters of the inhomogeneities, is sufficiently small, that these inhomogeneities are disjoint, 
and that their distance to R"* \ SI is larger than cq/2. Let /io and Eq denote the permeability 
and the permittivity of the background medium, and assume that /xq > and Sq > are 
positive constants. Let fij > and Sj > denote the permeability and the permittivity of the 
j-th inhomogeneity, Zj + aBj, these are also assumed to be positive constants. Introduce the 
piecewise-constant magnetic permeability 

1^ /ij, X G Zj + aBj, J = I . . .m. 

If we allow the degenerate case a = 0, then the function ^loi^) equals the constant fiQ. The 
piecewise constant electric permittivity, ea{x) is defined analogously. 

Consider solutions to the time-harmonic Maxwell's equations with exp(— iwt) time depen- 
dence. Let Ea be the electric field in the presence of the inhomogeneities. It solves the Helmholtz 
equation 

V • {—VEa) + uj'^EaEa = in n, 

with the boundary condition Ea — f on dfl, where w > is a given frequency. The electric 
field, Eq , in the absence of any inhomogeneities, satisfies the following equation: 

AEo + k^Eo = inn, (3) 

where fc^ = w^/igeo: with Eq = f on dfl. In order to insure well-posedness (also for the a- 
dependent case for a sufficiently small) we shall assume that k'^ is not an eigenvalue for the 
operator —A in L^{n,) with the Dirichlet boundary conditions. It has been shown in that 
the following asymptotic formula holds uniformly on dfl 

dEa, , dEo f dEa dE. , dG{x,y) 

av ov Jqq ov ov ov 

where the remainder o[a'^) is independent of the set of points {zj}JL^ provided that (|l|) holds, 
G{x, y) is a free space Green's function for A + k'^, and each Mj is a d x d, symmetric, positive 
definite matrix associated with the j-th inhomogeneity, called the polarizability tensor, which is 
given by 



\B,\Sa^+(^-l) [ yi^day, (5) 
VMo / JdB. ovj 



where, for 1 < I' < d, <j)i (y) is the unique function which satisfies 

A(/)r = in Bj and Tl'^\Ri, 

1 dcp^ ld^+ _ 1 

— I'j ■ ei' on O-Dj, 



with (l>i' continuous across dBj and Um|j^|^oo 4'i'{y) — 0- Here {ei'}f,^i is an orthonormal basis 
of R"*, Vj denotes the outward unit normal to dBj, superscripts — and + indicate the limiting 
values as the point approaches dBj from outside Bj , and from inside Bj , respectively. 

Our goal is to identify the locations {zj}JLi a-nd the polarizability tensors {MjYjLi of the 
small inhomogeneities Ba from the boundary measurements of on a given part Fi CC 
dfl. This inverse problem is more complicated from the mathematical point of view and more 
interesting in applications than the one solved in j^], because in many applications one cannot 
get measurements on the whole boundary. As in [Q, we want to reduce this reconstruction 
problem to the calculation of an inverse Fourier transform. Our work is a generalization of 
to the case of the reconstruction of small inhomogeneities from the measurements on a part of 
the boundary. 

In |lO| a different method is proposed for finding small inhomogeneities from the scattering 
data. In Q it is proved that the singular (delta-type) potentials are uniquely determined by 
the scattering data. In a numerical realization of the method proposed in | pO| is given. In 
Q a method for finding small inhomogeneities from tomographic data is given. In |^ analytic 
formulas are derived for electric and magnetic polarizability tensors for bodies of arbitrary 
shapes. In ||^ a number of multi-dimensional inverse scattering problems are studied and some 
of the methods developed there are related to our approach. 



2 Identification procedure 

Before describing our identification procedure, let us introduce the sets N{il) — {v : v ^ _ff^(f7)n 
H^{n),Av+k'^v = in 0} and N{n) = {v: v e H^{Ti)nH^{n), Av+Pv = in 17, w = on Fa}, 
where F2 = d^l \ Fi, where Fi is an open in dil subset. 

The general approach we use to recover the locations and the polarizability tensors of the 
small inhomogeneities is to integrate the solution Ea against special test functions in the set 
N{n). Let V be any function in N{n). As in ||l|, the following estimate can be derived from (^): 

/ I ds = a'y^{l-l^)VE,(z,).M,{^)Wv[z,) 

(6) 

in ^ ^ 

+ a^k^y^{l - ^-^)\B,\E^{z,)v{z,)^o{a^), 
3=1 

where \Bj\ stands for the volume of the set Bj. We want to make suitable choices for the test 
functions v in N{il) and the boundary condition Ea\dn in order to get simple equations for the 
unknown parameters, namely, for the points {zj}JLi and matrices {Mj}JLi. Similar idea was 
used and the associated numerical experiments have been successfully conducted in the case of 
the conductivity problem Q with boundary measurements on all of dfl. 

Let us describe our inversion method. Take r/ to be a vector in R"*, 77-*- a unit vector in 
R*^ which is orthogonal to i], and 7 a complex number. Then e'(''+''''' ^'^ is a solution to the 
Helmholtz equation in R'^ if and only if 7^ = A;^ — \ri\'^ and in this case e^^^~^^ ^'^ is also a 
solution to the Helmholtz equation in R''. For simplicity, let us consider the case where all the 
Bj are balls. In this case all the matrices Mj arc multiples of the identity matrix which makes 
our analysis simpler. 
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If is known on the whole boundary dQ then taking Ea = e^^'^'^'*'^ ^'^ on dVL we know 
from [y that 

^^"-(y)e'(''-^''")-^ ds{y)- I ^{e^n-iv^yv)E^{y) ds{y) 



an diy Jga dv 

j^i '-MM e J 

The main difficuhy in generalizing this approach to the case when is known only on a part 

Fi CC dfl is to construct a function Wa{x) in N{^1), that is asymptotically e*'''"'*'''^'''^ as a 
approaches 0. The following lemma holds. 

Lemma 2.1 Let Q,' CC be a -domain. Let rj S R'^ and t]-^ be a unit vector in R"^ that is 
orthogonal to rj. There exists Wa G N{Q) such that 

Wa{x) = e*(''-)")^)-^ + o(a'') and \/Wa{x) = 1(77 - jj^^je'^^-n^)-^ + o(a''), 

uniformly in fi'. 

This lemma is an immediate corollary of the following general density result. 
Proposition 2.1 The set N{^1) is dense, in the L^(f7') norm, in the set iV(il). 

Proof. Assume the contrary and let v G N{Q) be an element which cannot be approximated 
in L^(n') by the functions from N{U) with a prescribed accuracy. Then there is an element in 
N{n), which we denote again v, such that J^, vwdx = 0,\/ w £ N{n). Let Go be the Dirichlet 
Green's function in fi: 

AGo + fc^Go = Sy{x) in f^, 

Go = on an. 

Define i?3/2(-p^| = {p G H^/'^{Ti) such that there exists p G i/3/2(9f7),p|r2 = 0,p|ri = p}. 
Since any w e Af(ri) can be represented as follows 

/ ^^"^.""f^ P(y) dsiy),x&n, 

where p G -ff'^/^(ri) is arbitrary, we have 

f dGn{x,y) n w ^ 

Define := /j^, v{y)Ga{x,y) dy. Then = 0, |^|ri = 0, and 

A 1 2 f u in ri', 
^" + ^"=1 inf]\fr. 

By the uniqueness of the solution to the Cauchy problem for the Helmholtz equation, it follows 
that u = in \ Q', and so it = ^ = on dfl' . Since Au + k^u = w in fi', it follows by 
multiplying this equation by v and integrating by parts over fi' that J^, v'^dx = 0. Thus, v = 
in n', and, by the unique continuation, v = in which proves Proposition 2.1. In the above 
argument we assumed that u and v are real valued. This can be done without loss of generality 
since fc^ > 0. □ 
From Proposition 2.1 it follows that the function Wa in Lemma 2.1 can be constructed 
(numerically) by solving the minimization problem: 

min 11/ ^^^^p{y)ds{y)^e:(^-^^^y-U.^n'y (7) 
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We can take 



Wa{x) 



dGo{x,y) 



where Pa is any function in H^^'^iVi) such that 



Pa(y) ds{y), 



II I ds{y) - e^(''--''")-|U.(n') = o{a') 



(8) 



Since (A + k'^){wa — e*^'' -'■^) = in fi, by the standard cUiptic interior estimates 
obtain from (||) th at \\w a — e^(v-in ^'^H^Lfs^s^//) is of order o(a'*) for any s and any C^-( 
Jl" CC ri'. Lemma 2J follows then immediately from the Sobolev imbedding theorems. 

Now, if we choose Ea = e'^'^^'^'^ ^'^ on 9ri and v = Wa in then, since the points {zj} 
away from the boundary dfl, it follows from 
expansion holds: 



are 



and Lemma 2.1 that the following asymptotic 



Aa(?7) = 



dv 



Wa ds 



an 



dwa 



Ea ds 



^dj^e^'T'^ (1 



Mo 



Mo 



£0 



(9) 



Recall that e^*'''^^ (up to a multiplicative constant) is the Fourier transform of the Dirac function 
6-2z (a point mass located at —2zj). Multiplication by powers of rj in the Fourier space 
corresponds to differentiation of the Dirac function. Therefore, the function Aa{ri) is the inverse 
Fourier transform of a distribution supported at the points Zj. A numerical Fourier inversion 
of a sample of Aq,(?7) will yield the points Zj with a small error as a ^ 0. It is natural to use a 
fast Fourier transform for this inversion. One can estimate the number of the sampling points 
needed for an accurate discrete Fourier inversion using the Shannon's theorem j^, page 18. This 

number is of order (-r)^. One needs this amount of sampled values of rj to reconstruct, with 



resolution S, a collection of inhomogeneities that lie inside a square of side h. Once the points 
{^ilJLi are found, one may find {Mjyj^i by solving appropriate linear system arising from the 
asymptotic formula (^. If Bj are general domains, our calculations become more complicated, 
and eventually we have to deal with pseudo-differential operators (independent of the space 
variable x) applied to the same Dirac functions. In view of the asymptotic results derived in 
a similar approach may be applied to the full Maxwell equations in the presence of small 
dielectric inhomogeneities. 
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